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Abstract 

' It is shown that QCD with vector gluons predicts drastically different 

^| ■ multiplicity distributions compared to the scalar Ac^-model (in 6 dimen- 

sions). In particular, the QCD predicted minimum of cumulant moments 
' of the distributions at ranks about 5-^7 does not reveal itself in the scalar 

Ph ■ model. The oscillations at higher ranks survive only. Thus, the experimen- 

y—i . tal fact of that minimum at ranks about 5-^7 supports the vector nature 

of gluons and singular kernels in QCD and gluodynamics. 

> : 

in 

O 1 Introduction 

<n : 

The multiplicity distribution of particles (partons) in QCD (or gluodynamics) 
has been a long-standing problem. The solutions of equations for generating 
functions of multiplicity distributions have been obtained in the lowest double 
logarithmic approximation |]T| . The shape of the distribution appears to be com- 
pletely different from that observed in multiparticle production experiments for 
secondary hadrons, even though the energy dependence is rather reasonable and, 
moreover, the shape fulfils the so-called KNO-scaling ||. The theoretical shape 
is much wider. In terms of the moments of the distribution it implies much faster 
^ ■ increase of them with rank increasing, compared to experimental data. 

Meantime, the theoretical attempts to take into account the conservation laws 
more precisely were published 0, f|, [|. The multiplicity distribution became 
rather narrow||. However, in view of some additional assumptions it is difficult 
to estimate how rigorous the results are. Moreover, the theoretical moments of 
the distributions differ from experimental ones. 

The problem was solved when the solution of the equations for generating 
functions was found at higher (than double-logarithmic one) approximations in 
gluodynamics [|, |7| and QCD || in case of the running coupling constant. Later, 
the exact solution of these equations at fixed coupling constant was obtaned|J. 
The most remarkable prediction M of these findings is the minimum of the cumu- 
lant moments of the multiplicity distribution at the values of ranks close to 5 and 
subsequent oscillations of those moments at higher values of the rank [[7| which 
have been confirmed in experiment at analysis of hadron multiplicity distribu- 
tions in multiparticle production at high energies for various colliding particles 
[ |I0|j (see also the review paper [|TT|j). 
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Let us note also that no phenomenological distributions well known and widely 
used in physics (i.e. those distributions of probability theory as Poisson, nega- 
tive binomial, geometric, fixed multiplicity) can reproduce such a behaviour. In 
particular, the negative binomial distribution gives rise to the positive and mono- 
tonically decreasing ratio of cumulant to factorial moments, while QCD predicts 
the negative values of this ratio at the first minimum and its subsequent oscilla- 
tions about the abscissa (where the ranks of the moments are plotted). Even the 
so-called modified negative binomial distribution, specially designed empirically 
12] to get the best fit of experimental multiplicity distributions has failed in some 
p~3]| to produce the relevant minima and oscillations of the moments. Thus 
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the analysis of multiplicity distributions in terms of cumulant and factorial mo- 
ments as well as of their ratio is at present the most sensitive method of revealing 
the specific typical details of these distributions. 

The very existence of the minimum in QCD and its location are closely defined 
by the singularity of the gluon kernel (the vector nature of gluons) and by the 
value of the coupling constant in QCD. In that connection it is of interest to 
understand how important the both factors are. According to the formulae of the 
paper || , the main contribution to the minimum location is inverse proportional 
to the QCD anomalous dimension (i.e. to the coupling constant) and therefore 
varies in QCD in the limits of its "running" property i.e. comparatively slow 
and hard to observe if, especially, one takes into account that the ranks of the 
moments can be integer only. The more important property is the vector nature of 
massless gluons giving rise to the singularity of the kernel of the equation for the 
generating function. That is why it would be desirable to confront the QCD (or 
gluodynamics) predictions to those of the scalar fields model. Fortunately, there 
is such a model [[b|] possessing, besides others, the property of the asymptotical 
freedom what helps to reduce the impact (even if it is not very essential) of the 
problem of the value of the coupling constant. This is the A<^ 3 -model in the 
6-dimensional space-time. 

Therefore the purpose of the present paper is to get the knowledge of the 
behaviour of the multiplisity distribution moments for the A(^ 3 -model in the 6- 
dimensional space-time and to compare the obtained results with those of QCD 
and of the simplest phenomenological distributions of the probability theory. The 
main conclusion which we have got from it is that the vector nature of gluons 
is very crucial for the shape of the multiplicity distribution and qualitatively 
changes its moments behaviour so that QCD confirms its predictive power once 
again even at purely partonic level while the scalar model has been unable to 
reproduce the qualitative features of experimental data. 
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2 Equations and their solutions 



The theoretical problem, we are addressing at, is as follows. Let a scalar, strongly 
virtual particle (with high time-like 4-momentum squared) has been produced in 
some collision. During its evolution the virtuality decreases due to emission of 
other scalar particles. What is the multiplicity distribution of secondary par- 
ticles if the emission is controled by the A<^ 3 -interaction in the 6-dimensional 
space-time? This problem is analogous to the production and evolution of the 
pair of strongly virtual quark and antiquark in e + e~-annihilation. The increased 
dimensionality of the space-time, as we have mentioned, is necessary here just 
to get the asymptotical freedom in that theory and to make it closer to the real 
situation in QCD (or in gluodynamics). 

The evolution of such a "scalar jet" is described by the equation for generat- 
ing functions (or functionals) analogous to the commonly used "birth and death 
equation" and different from QCD equations by its kernel, describing the inter- 
action vertex in the theory. However, just this distinction strongly influences the 
multiplicity distribution and its moments. The principal difference is that the 
QCD kernel is singular while in the scalar model it is regular. Physics corollary 
is an approximate equipartition of the "parent" energy among its "children" in 
the scalar model while in QCD the energy is shared in unequal parts i.e. one of 
the produced partons used to get much higher energy than another one. 

Let us turn to our problem and remine some general notations. The generating 
function is 

oo 

G(y,z) = £(l + *) B P„(l/), (1) 

n=0 

where P n (y) is the multiplicity distribution, y = hn(Q/Q ), Q is the jet virtuality, 
Qo =const. 

The normalized factorial moments (F q ) and cumulants (K q ) are defined by 
the generating function as 

oo q oo q 

G{z) = £ -<n)<% \nG(z) = £ Z -(n)«K q , (2) 

q=0 1- q=0 

where (n) is the average multiplicity of particles (partons). 

The generating function satisfies the non-linear integro-differential equation 
[|TJ (the prime denotes the ^/-derivative) : 

G'(y) = f 1 d^K(Ol 2 [G(y + \n£)G(y + ln(l - 0) " G(y)), (3) 
Jo 

where K(£) is the kernel of the equation and in our cases is written as: 
1) for gluodynamics ]!]: 

a-(0 = i/£-(i-0[2-£(i-0L ( 4 ) 
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2) for A<^ 3 -model in 6 dimensions ||T4]| : 



^(0 = 6^(1-0- (5) 

7o = 2N c a s /n, a s = 2ir/lly in QCD. In Ay? 3 -theory one does not have any 
physics normalization and has to rely on QCD 7q in choosing the numerical 
value of A. 

Now we discuss Ay2 3 -model directly and consider first the running coupling 
case, following the method of the approximate solution used in |J,@ for gluody- 
namics. After the Taylor series expansion of the generating function at point y 
has been done, one gets 



G\V) = ll{G{y){G{y)-l)+G{y)Y.^) n hnG^\y)+ £ {-l) n+m h nm G^\y)G^\y)}, 

n=l n,m=l 

(6) 

whereQ 

h n = ^(-1)" H dm - In»(l -0 = 12(^r - 3^), 
l) n+m r^l-O ln"(0ln m (l-0 

J 
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n\m\ Jo 
or within the approximation used in 0, one gets 



(lnG(y))' = 7o 2 {G(y) - 1 - hG'iy) + h 2 G"(y) + h xl G'{y){\nG(y))'}, (7) 

where h x = 5/3, h 2 = 19/18, h n = (37 - 3tt 2 )/18. 

The anomalous dimention 7 of the A<£> 3 -model is defined similar to QCD : 
(n) = exp(J y j(y')dy') and we substitute (||) to (0). The coefficients in front 
of z^-terms should be equal on both sides wherefrom the moments satisfy the 
equation 

^ = t^ho g ^ |T + v[ T + 7')' (8) 

x = g 7 , = 7o 2 [l/^ - 5/3 + /i 27 7 7 ], k q = K q /{q - 1)!, f q = F q /(q - 1)! and we 
have used the relation 

f q = k q + q ^^. (9) 
1=1 1 



1 We are using the same notations hi for the coefficients in the Ay 3 -model as in gluodynamics 
: (t|, and hope that it does not produce any misunderstanding even though their numerical 
values differ in the both cases. 
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The recurrent formula (p|) has been used to calculate k q and to get f q with the 
help of eq.@ as well as their ratio H q if 7 is known. 

The equation for 7 obtained from @ in the similar way to (0) but for q = 1 
is written as 

7-7o[l-^7 + ^(7' + 7 2 )] = 0, (10) 
what represents 7 in terms of 70 as 

7 = 7o 2 [l-^7o 2 ]+0(7o 6 )- 

The results of computing H q according to (|8|) for 70 = 0.48 are shown in 
Fig.l. In distinction to QCD results obtained within the same approximation 
in 57, there are still no oscillations and H q tends to constant asymptotics (it 
reminds of QCD in the lower order in qj - see f|). 

Varying 70 within the wide limits, one does not observe any essential quali- 
tative changes. The minimum is slightly shifted to larger values of q at smaller 
7o- 

Thus we notice that the absence of the singular term in the kernel @ in 
the A(^ 3 -model compared to the gluodynamics (the formula (^)) gives rise to 
qualitatively different behaviour of the moments. 

In the case of the fixed coupling constant, the equation (|3|) can be solved ex- 
actly . Assuming that the ^/-dependence is contained completely in the average 
multiplicity (n) and is of the kind (n) = exp^y), one gets the recurrent formula 
for the factorial moments 

6 7o - 
fq = x{2 + x)(3 + x) + 7o 2 (x - l)(x + 6) § fq ' lfh (U) 

where f q = f q T(x + 2)/q = F q T(x + 2)/q\ and 

zj 1 H q _[ fn—lfl 

tiq=l- 2^ 



1=1 
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The analogous formula in gluodynamics looks like (it can be easily obtained 
from the formulae of || if the quark degrees of freedom are neglected) 

7o 9-1 

fq = x-'lM^^- 1 ^ (12) 

where 

M g = i + *(l)-*(l + x) + ^- 2 ' 



x w v ' 12 1 + x (2 + x)(3 + x) 

( 1+x \ 1 1 

*' ~ I ~h~ J (l + h)(l + (q-lh) + 2(2 + x)(3 + x) 
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In Figs. 2 and 3 we show the results of computing according to formulae (O) 
and ( pp correspondingly. One can conclude from them that the influence of the 
singularity in the kernel K(£) on the location of the minimum consists in its 
shift to smaller values of ranks q ( from q\ w 14 to q2 ~ gi/2 = 7) with the 
accompanying diminishing of the "period" of the oscillations. 



3 Conclusion 

The main conclusion one gets from above consideration is that the scalar \ip 3 - 
model in 6 dimensions with approximate equipartition of energy among the pro- 
duced particles gives rise to the qualitatively different predictions about the be- 
haviour of the moments of multiplicity distributions compared to the gluody- 
namics or/and QCD, where the vector nature of massless gluons implies the 
singularity of the kernel of the equation and, therefore, the drastically unequal 
shares of energy for produced particles. This fact provides the minimum in the 
ratio of cumulants to factorial moments at ranks near 5 -j- 7 in gluodynamics 
and QCD, while in the scalar model such a minimum does not appear and the 
oscillations at higher values of ranks survive only. 

From the purely theoretical approach, this conclusion is hardly extremely im- 
portant. However, it becomes essential when one confronts it to the experimen- 
tally known facts about hadron multiplicity distributions in real events [K], |TT|] . 



It happened that the minimum of H q at q ~ 5 was observed [1C] in experimental 
data just after its prediction |J in gluodynamics. Even though this prediction 
was done about partons and not about final hadrons, the very fact of the presence 
of the minimum with the same location as predicted looks very impressive. In 
combination with the inability of the phenomenological models as well as of the 
theory field model considered above to reproduce the existence of the minimum 
and its location, it confirms once again our belief in the predictive power of QCD. 
Other numerous facts (the hump-backed plateau of the rapidity distribution, the 
heavy quark effects etc. - see in detail in U and in later papers) support the 
conclusion that the qualitative effects of QCD predicted at the parton level find 
out their confirmation by the experimental data, and Monte-Carlo simulations 
provide the quantitative estimates of the influence of hadronization. 
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